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ABSTRACT 

We explore the possibility of detecting primordial non-Gaussianity of the local type using weak 
lensing peak counts. We measure the peak abundance in sets of simulated weak lensing maps corre- 
sponding to three models /nl = 0, —100, 100. Using survey specifications similar to those of EUCLID 
and without assuming any knowledge of the lens and source redshifts, we find the peak functions of 
the non-Gaussian models with /nl = ±100 to differ by up to 15% from the Gaussian peak function 
at the high-mass end. For the assumed survey parameters, the probability of fitting an /nl = peak 
function to the /nl = ±100 peak functions is less than 0.1%. Assuming the other cosmological param- 
eters known, /nl can be measured with an error A/nl ~ 13. It is therefore possible that future weak 
lensing surveys like EUCLID and LSST may detect primordial non-Gaussianity from the abundance of 
peak counts, and provide complementary information to that obtained from the cosmic microwave 
background. 



1. INTRODUCTION 

The inflationary paradigm is the leading theory of the 
early Universe, of fundamental interest for cosmology 
and particle physics. Understanding the mechanism and 
energy scale of inflation remain major goals to attain, de- 
spite the continuous and fervent efforts invested in this 
field. 

A measurement of primordial gravitational waves 
would pin down the energy scale of inflation, though it 
still belongs to the not-so-near future. One possible way 
to discriminate between single- and multi-field inflation 
models is to test the Gaussianity of the pr i mordi al den- 
sity fluctuations Creminelli & Zaldarriaga (2004). The 
Cosmic Microwave Background (CMB) has been so far 
the main and cleanest inflationary probe. Recent re- 
sults fro m the Wilkinson Micro wave Anisotropy Probe 
(WMAP) IKomatsu et~aTl (|2010h suggested tentative la- 
level evidence of primordial non-Gaussianity of the local 
type, defined by the equation: 



*(x) = <j>(x) + /nl[0 2 (z) - {<f{x))]. 



(1) 



The parameter /nl quantifies the local quadratic devia- 
tion of the Bardeen potential $ from a Gaussian poten- 
tial 4>, and it is c urren tly constrained to the value 32 ± 21 
IKomatsu et~aL1 (|2010D . 

It has long been su g gested iGrinstein fc W ise 

iFrv fc Scherrerl 



Lucchin & Matarrcse (1 



1986) 



1994); 



Matarrese et alj (|2000D that low-redshift observables can 
also be used to measure primordial non-Gaussianity, de- 
spite the fact that the density field at such redshifts 
is strongly non-Gaussian due to the action of grav- 
ity. In the local non-Gaussianity models, there are 
mainly two effects on low-redshift observ ables as re- 
cently outlined in the comprehe nsive study oflDalal et alj 
(120081). and further explored in Desj acques et al.|( |2009); 
IGrossi et all (120091) ;JptQepich et alJ(j2010l) ; ILam fe Shethl 
(|2009D : ISmith et all (|2010D . First, /nl induces a scale- 



dependence in the bias of dark matter halos, which af- 
fects primarily the largest scales, i.e. k < 0.02/iMpc -1 . 
Thus one can in principle separate the /nl scale de- 
pendence from the gravitational one, which occurs on 
smaller scales. Second, the abundance of massive ha- 
los is higher /low er for positive / negativ e values of /nl 
iMatarrese et all (|2000fk ILo Verde et all (|2008D . As a fi- 
nal note, we emphasize the importance of high- and low- 
redshift constraints on non-Gaussianity as their compari- 
son may prov i de ins i ght into the scale-depen dence of /nl 
IBvrnes et all (|2010D ; [Shandera et alj (|2010l) . 

In this paper we shall numerically investigate the sen- 
sitivity of weak gravitational lensing (WL) peak counts 
to primordial non-Gaussianity of the local type. The po- 
tential of WL surveys to con s train / nl has alrea d y bee n 
tackled lAmara fc Refregierl (120041): IPace et all 
iFedeli fc Moscardinil (|2010l ); IQguri fc Takada 
though without considering shear peaks. Peak counts 
are a natural candidate for /nl studies, since the 
largest of them are caused primarily by massive ha- 
los. For W L studies of peak st a tistics, we mention 
the w o rks of IHamana et aLl (120041); |Hennawi fc Spergel 
(l2005l);IDe Putter fc WhitT(l2005l);lMarian et alj (12009) ; 



ap (120101) : 
M 12010), 



Kratochvil et aLl (120101) ; iDietrich fc Hartlad (|2O10l) : 
Marian et al.1 (|20"lC ) and the references therein. Should 
peak counts prove to be a sensitive /nl probe, then 
could easily use large future surveys like EUCLID 



one 
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easi ly 

IRefregier et al.l (120101) or the Large Synoptic Survey Tele - 
scope (lsst) [LSST Sci ence Collaborations et all (|2009T ) 
to obtain low-redshift constraints on primordial non- 
Gaussianity. 

2. METHOD 

The observable that we use is the convergence field, 
i.e. the matter density projected along the line of sight 
and scaled by a geometrical factor. We study simulated 
WL convergence maps created from ray-tracing through 
a suite of A-body simulati ons, generated w ith the pub- 
licly available code gadget ISprin gcl ( 200l). A subset of 
these simulations was used and described in the work of 



2 



iDesiacaues et a l. ( 2009) . Three values of /nl are consid- 
ered: 0, -100, +100, while the other cosmological pa- 
rameters are ke pt fixed. The cosmo logy matches the 
WMAP5 results Komatsu et all (|2009f ). We have a total 
of 18 simulations, with 6 realizations per /nl value. The 
initial conditions for each of the 6 sets of /NL-model re- 
alizations are matched to reduce the cosmic variance on 
the comparison of the peak functions corresponding to 
each model. The box size is 1600 h^ 1 Mpc, the number 
of particles is N — 1024 3 , and the softening length is 
'soft = 40/i _1 kpc. 



We consider a surv ey similar to EUCLID | Rcfrcgic r et al 
(|2010[) and to lsst ILSST Science Collaborations et al 



(2009) for the WL simulations with: an rms <r 7 = 0.3 for 
the intrinsic image ellipticity, a source number density 
"■gal = 40arcmin -2 , and a redshift distribution of source 
galaxies given by V(z) = 1.5 z 2 /zq exp[— (z/zq) 1 - 5 ], 
where zq — 0.6. The median redshift of this distribu- 
tion is z mo d = 0.86. 

From each iV-body simulation we generate 16 indepen- 
dent fields of view. Each field has an area of 12 x 12 deg 2 
and is tiled by 4096 2 pixels, yielding an angular resolu- 
tion 8 pix — lOarcsec and a total area of ~ 14000 deg 2 
for each /nl model. The effective convergence k in each 
pixel is calculated by tracing a light ray back through the 
simul ation with a Multiple - Lens- Plane ray-tracing algo- 
rithm |MbirtetaI| ([2002H IHl). Gaussian shape noise 
with variance cr 2 /(n ga i # 2 ix ) is then added to each pixel, 
which creates a realistic noise level and correla tion in the 
filtered convergence field IHilbert et~aTl (12007% 

For the peak finding we use an aperture filter Schneider 
(1 19961). matching an Na varro-Frenk- White (NFW) pro- 
file lNavarro et all (|1997| ) convolved with a Gaussian func- 
tion. Thus we adopt the convergence model: ft mo dci = 
k nfw o /causs, where /Gauss is a Gaussian function of 
width / x i so ft- 'soft is the softening length of the sim- 
ulations, / = 1.5 for M < 7 x 10 14 /i _1 M s and / = 2 
otherwise. Here the NFW conve rgence profile 

trunca ted at the virial radius, defined in iMarian et al.l 

(2010) . This model agrees very well with the measured 
convergence profiles of the peaks in the maps. It is use- 
ful both when working with simulations, since it accounts 
for the lack of resolution below the softening scale, and 
also when using real data, since shear data is difficult to 
obtain near the centres of clusters. 

The amplitude of the smoothed field at a point Xq is 
given by: 



M(x ) 



d 2 x W(x — x)k(x), 



(2) 



where W is our filter and k is the convergence field. The 
filter W can be written as follows: 



W(x) = C v 



K model(2') K model(-R) 



o- 2 /n 



gal 



(3) 



where C w is a normalization constant and R is the aper- 
ture radius, i.e. the radius over which the filter is com- 
pensated, and k(x) = 2/x 2 J Q dyyK(y). We choose the 
normalization constant to be: 



7 
n S al 



M N 



Jdh 



model 



(x) 



(4) 



If Xq is the location of a peak created by an NFW cluster 
of mass M NFW and redshift z, and the convergence field is 
smoothed with a filter tuned to precisely such a cluster, 
i.e. Kmodci in Eq. ((3]) corresponds to the same M NFW and 
z, then this filter returns a maximum S/N at Xo. At this 
location, the amplitude of the smoothed map is: 



M(x \M NFyv ,z) = M N 



(5) 



The peak is assigned the mass M NFW . If the peak does 
not correspond to a real halo, and it is the result of line- 
of-sight projections, then it can still be assigned an 'ef- 
fective' mass. In practice, we smooth the convergence 
field with filters of various masses, which yield different 
amplitudes (larger and smaller than the filter mass) at 
the location of peaks. We interpolate these amplitudes 
to determine the filter mass that would satisfy Eq. ([5]). 

We choose R to be the virial radius of the cluster to 
which the filter is tuned: R = (3M NFW /800 np) 1 ^ 3 , where 
p is the mean density of the universe. Thus a filter at a 
given redshift can be specified either through the mass 
M NFW or through its siz e R. We adopt the mas s con- 
vention of Sheth-Tormen ISheth fc Tormei] (fl999h . with 
an overdensity defined as 200 x the mean density (not 
the critica l densi ty) , and the concentration parameter of 
iGao et al.l ()2008[ ). We evaluate the S/N of the mass esti- 
mator in a simplified scenario where we ignore projection 
effects, and consider only the intrinsic ellipticity noise. 
In this case the variance of the estimator is given by: 
Var(M) = M NFW C W , and the S/N is obtained by com- 
bining this with the above equations. Note that while 
the mass estimator does not depend on the shape noise 
due to the normalization constant C, the S/N scales with 
it as aperture filters usually do: S/N ~ (n ga i/ V 2 ) 1 / 2 . 

Finally, in the case where no shape noise i s included, we 
use th e same filter as in our previous work IMarian et aLl 
<|2010f >. This filter can also be obtained by formally tak- 
ing cr 2 /n ga i = 1 in Eqs (©, Q above. 

The analysis of the convergence maps is carried out in 
two situations: with and without shape noise. We adopt 
a very conservative approach in which we consider known 
only the redshift distribution of the source galaxies and 
the shape noise level, without any other information on 
the sources or on the detected peaks. We also do not 
resort to tomographic techniques. While this is an overly- 
pessimistic scenario for a next-generation lensing survey 
like EUCLID or LSST, our goal here is to provide a proof 
of concept of the possibility of using WL peak counts 
to constrain primordial non-Gaussianity, rather than the 
final quantitative answer to this question. 

We perform a hierarchical smoothing of the maps with 
filters of various sizes, from the largest down to th e small- 
est, as described in IMarian et all (|2009L 120101 ) . This 
approach removes the problem of 'peaks-in-peaks' and 
it also naturally eliminates the dependence of the mea- 
sured peak function on a particular filter scale. Since 
the median redshift of the source distribution is 0.9, we 
adopt a fixed redshift of 0.3 for the matched filter de- 
scribed above. For this redshift, the scale of the fil- 
ter varies from corresponding masses of 2 x 10 15 ft. -1 M Q 
to 3 x 10 13 h^ 1 M Q in the absence of shape noise, and 
10 14 hr 1 Mq in the presence of it. The latter lower-limit 
choice of the filter scale is due to the fact that shape 
noise contaminates seriously the smaller peaks; imposing 
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Fig. 1. — Upper panels: The measured peak functions for /nl = +100 (red circles) and /nl = —100 (blue triangles) for a 
survey with a median redshift z mc( j = 0.9 and an area of 14000 deg 2 . The filter used corresponds to redshift 0.3. The peak 
functions in the left panel are measured in the absence of shape noise, while in the right panel shape noise is included (as- 
suming 40 galaxies/arcmin 2 ), and only the peaks with a S/N > 3 were selected (hence the smaller mass range on the ic-axis). 



Lower panels: The fractional difference of the peak functions for the /nl = ±100 models and the Gaussian model (red circles/blue 
triangles), without/with shape noise (left/right panel). The points are obtained as an average over all the fields of the fractional difference 
of the /nl i 100 peak functions in every field and the average peak function in the Gaussian model. 



a minimum S/N threshold alleviates but does not remove 
the contamination. Peaks are detected in the smoothed 
maps as local maxima, and are assigned a mass as de- 
scribed above. 

3. RESULTS 

Figure [T] presents the main result of this work. The 
upper panels illustrate the peak functions measured in 
the /nl = ±100 cosmologies, in the absence (left panel) 
and presence of shape noise (right panel) . For each /nl 
model, the points are the average of the peak counts mea- 
sured in all 96 fields and the error bars represent errors 
on the mean. In the case of shape noise, we select only 
peaks with a S/N > 3. We note that the two average 
peak functions are clearly distinct at the high-mass end, 
i.e. for M > 3 X 10 14 /i _1 M Q . This is better seen in 
the lower panels which show the difference of the peak 
abundance measured in the /nl = ±100 cosmologies, 
relative to the /nl = peak abundance. We have min- 
imized the impact of the matched initial conditions of 
the simulations on the error bars: the fractional differ- 
ence is computed as an average of the /nl ± 100 peak 
abundance in each field ratioed to the average Gaussian 
peak abundance of 96 fields (which represents our best 
model for the true Gaussian peak abundance), as op- 
posed to the average of the ratio of the non-Gaussian 
and Gaussian peak functions measured in the same field. 



The latter would have removed the cosmic variance of 
the fractional difference, because of the matched initial 
conditions of the simulations from which we built the 
convergence maps. 

Just like in th e case of the 3D halo mass function, 
see f or example iDesjacques et al.1 (|2009l ): iSmith et al.l 
2010)), the peak functions for the /nl models show a 
deviation from the Gaussian case. Unlike the 3D studies 
which have presented halo mass functions measured at 
a single redshift, the peak functions that we show here 
combine peaks in the redshift range of the source dis- 
tribution, and therefore are not as regular and symmet- 
ric as their 3D counterparts. The asymmetry is most 
likely due to modifications of the S/N of peaks by line- 
of-sight projections, and also by shape noise contami- 
nation. The trend is similar however, with high-mass 
peaks displaying the largest deviation. For both /nl 
models, this is about 10-15% for the largest mass bins 
i.e. M > 4 x 10 u h- 1 M e . 

To quantify the significance of the deviation, we per- 
form a x 2 — test. For the fiducial model /nl = 0, we 
estimate the covariance of the counts in a field. We use 
the covariance of the mean to obtain the x 2 ■ For both 
/nl±100 we find a probability < 0.1% to fit the /nl±100 
peak functions with an /nl = peak function. This is 
also true if we consider only the diagonal elements (the 
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variance of the mass bins) instead of the full covariance 
matrix, and also if we vary the mass bins. We also use 
the measured counts to estimate the Fisher error that a 
14000 deg 2 WL survey would yield on /nl- Assuming 
all other cosmological parameters known, the forecasted 
error is A/nl ~ 13 for the fiducial value /nl = 0. In the 
above we considered that the peak abundance scales lin- 
early with /nl, similar to the 3D halo abundance. The 
values of the x 2 and the Fisher error are largely main- 
tained also if we minimize the impact of the matched ini- 
tial conditions of the simulations, by using the first three 
simulations to compute the /nl = +100 peak function 
and the last three for the /nl = —100 function. 

Though these results are already very encouraging, it 
is possible to improve measurements of primordial non- 
Gaussianity from WL surveys even if one considers only 
peak counts. The most important is probably the use 
of tomography, as numerical and analytical studies of 
the 3D mass functions have shown that deviations of 
the halo abundance from the Gaussian case significantly 
increase with redshift. Tomography would allow peaks 
to be separated not only in terms of their S/N-mass, 
but also of their redshifts, thus acquiring more sensitiv- 
ity to /nl- On the other hand, one has to beware a 
possible degeneracy with the amplitude of the matter 
power spectrum, er 8 . This can be solved by either using 
as a prior very g ood knowledge of as from other p robes, 
such as the CMB lThe Planck Collaboration! (|2006f ). WL, 
and large-scale structures or by combining several ob- 



servables sensitive to b oth /nl and as, as exemplified in 
lOguri fc Takadal (|2010l ) . Our immediate goals are to fur- 
ther study how WL can be used to constrain primordial 
non-Gaussianity, to build improved /nl estimators from 
WL observables, and to forecast /nl constraints based 
on these observables. 

For now we convey a simple, yet powerful state- 
ment: future WL surveys could detect primordial non- 
Gaussianity of the local type from at least one statistic- 
peak counts. In particular, surveys like EUCLID and LSST 
should be able to provide complementary information to 
the /nl infor mation obtained from the ongoin g CMB 
mission Planck lThe Planck Collaboration! (|2006D . 
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